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Abstract 

A strictly linear evolution of the cosmological scale factor is surprisingly 
an excellent fit to a host of cosmological observations. Any model that can 
support such a coasting presents itself as a falsifiable model as far as classical 
cosmological tests are concerned. Such evolution is known to be comfortably 
concordant with the Hubble diagram as deduced from data of recent super- 
novae la and high redshift objects, it passes constraints arising from the age 
and gravitational lensing statistics and clears basic constraints on nucleosyn- 
thesis. Such an evolution exhibits distinguishable and verifiable features for 
the recombination era. This article discusses the concordance of such an evo- 
lution in relation to minimal requirements for large scale structure formation 
and cosmic microwave background anisotropy along with the overall viability 
of such models. 

While these results should be of interest for a host of alternative gravity 
models that support a linear coasting, we conjecture that a linear evolu- 
tion would emerge naturally either from General Relativity itself or from a 
General Relativistic theory of a non-minimally coupled scalar field theory. 
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1 Introduction 



Large scale homogeneity and isotropy observed in the universe implies the 
Friedman-Robertson- Walker (FRW) metric: 

dv^ 

ds 2 = dt 2 - a(t) 2 \ —- - + r 2 {d0 2 + sin 2 9d(t) 2 )} (1) 
1 — Kr 

Here K = ±1,0 is the curvature constant. In standard "big-bang" cosmol- 
ogy, the scale factor a(t) is determined by the equation of state of matter 
and Einstein's equations. The scale factor, in turn, determines the response 
of a chosen model to cosmological observations. Four decades ago, the main 
"classical" cosmological tests were (1) The galaxy number count as a func- 
tion of redshift; (2) The angular diameter of "standard" objects (galaxies) 
as a function of redshift; and finally (3) The apparent luminosity of a "stan- 
dard candle" as a function of redshift. Over the last two decades, other tests 
that standard cosmology has been subjected to are: the early universe nucle- 
osynthesis constraints, estimates of age of the universe in comparison to ages 
of old objects, statistics of gravitational lensing, constraints from large scale 
structure formation and finally, the physics of recombination as deduced from 
cosmic microwave background anisotropy. 

In this article we explore concordance of the above observations with a 
FRW cosmology in which the scale factor evolves linearly with time: a(t) oc t, 
right from the creation event itself. The motivation for such an endeavor has 
been discussed at length in a series of earlier articles [1]. First of all, the use 
of Einstein's equations to describe cosmology has never ever been justified. 
The averaging problem and the continuum limit in General Relativity have 
not been properly addressed. In any case, most treatments have not consid- 
ered the retarded effects in their full generality[2, 3, 4]. Averaging Einsteins 
equations across horizon lengths lacks self consistency. On the other hand, 
Newtonian cosmology, applied to an exploding Milne ball in a flat space-time 
[see eg. [5, 6]] gives a unique linear coasting cosmology viz. the FRW [Milne] 
metric with a(t) = t. Such a cosmology does not suffer from the horizon 
problem. As a matter of fact, a linearly evolving model is the only power law 
model that has neither a particle horizon nor a cosmological event horizon. 
Linear evolution is also purged of the flatness or the fine tuning problem 
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[1, 7, 9]. The scale factor in such theories does not constrain the matter den- 
sity parameter. The Linear coasting characteristic of a Newtonian cosmology 
can be dynamically generated. As a matter of fact, it is a generic feature in a 
class of models that attempt to dynamically solve the cosmological constant 
problem [8, 7, 9, 10, 11]. Particularly appealing in a fourth order conformally 
invariant model of Mannheim and Kazanas [11]. Gravitation is determined 
by a fourth order (Weyl tensor squared) action that identically vanishes for 
the conformally flat FRW metric. A non- minimally coupled scalar field then 
produces an effective repulsive gravitation that quickly constrains the uni- 
verse to a linear coasting. Non-conformally flat perturbations around this 
linear coasting background can be effectively described by eqn(3) below. 

We may take any of the above as the basis for our linear coasting con- 
jecture. In what follows, we assume that an homogeneous background FRW 
universe is born and evolves as a Milne Universe about which a matter dis- 
tribution and standard General Relativity would determine the growth of 
perturbations. The following sections describe standard aspects of cosmol- 
ogy that follows. 



2 A linearly coasting cosmology 
2.1 Classical Cosmology tests 

Kolb[12] was probably first to demonstrate that data on Galaxy number 
counts as a function of red-shift as well as data on angular diameter distance 
as a function of red-shift do not rule out a linearly coasting cosmology. Un- 
fortunately, these two tests are marred by effects such as galaxy mergers and 
galactic evolution. For these reasons these tests have fallen into disfavour as 
reliable indicators of a viable model. 

The variation of apparent luminosity of a "standard candle" as a function 
of red-shift is referred to as the Hubble test. With the discovery of Supernovae 
type la [SNe la] as reliable standard candles, the status of Hubble test has 
been elevated to that of a precision measurement. Recent measurements by 
the supernovae cosmology project [13] eliminated the "minimal inflationary" 
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prediction defined by f2 A = and VLm = 1- The data has been used to assess 
a "non- minimal inflationary cosmology" defined by Q\ ^ 0, f2 A + = 1. 
The maximum likelihood analysis following from such a study has yielded 
the values tt M = 0.28 ± 0.1 and fi A = 0.72 ± 0.1 [13, 14, 15]. 

For a general power law cosmology with the scale factor a(t) = kt a , 
with k, a arbitrary constants, it is straightforward to discover the following 
relation between the apparent magnitude m(z), the absolute magnitude M 
and the red-shift z of an object: 

m (z) = M + 5logH + 5log(^r(l + z)kS[— * (^-(l - (1 + ^)] 

ri yl — a)K ti 

(2) 

Here S[X] = X, sin(X) or sinh(X) for K — 0, ±1 respectively, and M = 
M - 5log(H Q ) + 25. The best fit turns out to be a = 1.001 ± .0043, K = -1. 
[1]. The minimum x 2 per degree of freedom turns out to be 1.18. This is 
comparable to the corresponding value 1.17 reported by Perlmutter et al for 
the constrained non-minimal inflationary cosmology. Linear coasting is as 
accommodating for more recent high red-shift objects as the standard non- 
minimal inflationary model. The concordance of linear coasting with SNela 
data finds a passing mention in the analysis of Perlmutter [13] who noted 
that the curve for = Qm = (for which the scale factor would have a 
linear evolution) is "practically identical to bestfit plot for an unconstrained 
cosmology" . We display this "practical coincidence" in figure I that includes 
the more recent high redshift objects. 

The age estimate of the (a(t) oc t) universe, deduced from a measurement 
of the Hubble parameter, is given by t = (i? ) _1 . The low red-shift SNela 
data [16, 17] gives the best value of 65 km sec" 1 Mpc" 1 for the Hubble pa- 
rameter. The age of the universe turns out to be 15 x 10 9 years. This is ~ 50% 
greater than the age inferred from the same measurement in standard (cold) 
dark matter dominated cosmology (without the cosmological constant). Such 
an age estimate is comfortably concordant with age estimates of old clusters. 

A study of consistency of linear coasting with gravitational lensing statis- 
tics has recently been reported [1]. The expected frequency of multiple image 
lensing events probes the viability of a given cosmology. A sample of 867 high 
luminosity optical quasars projected in a power law FRW cosmology gives an 
expected number of five lensed quasars for a linear coasting cosmology. This 
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indeed matches observations. Thus a strictly linear evolution of the scale 
factor is comfortably concordant with gravitational lensing statistics. 

2.2 "The precision" tests 

a) The Nucleosynthesis Constraint What makes linear coasting particu- 
larly appealing is a recent demonstration of primordial nucleosynthesis not to 
be an impediment for a linear coasting cosmology [1, 18]. A linear evolution 
of the scale factor may be expected to radically effect nucleosynthesis in the 
early universe. Surprisingly, the following scenario goes through. 

Energy conservation, in a period where the baryon entropy ratio does not 
change, enables the distribution of photons to be described by an effective 
temperature T that scales as a(t)T = constant. With the age of the universe 
estimated from the Hubble parameter being « 1.5 x 10 10 years, and To ~ 
2.7K, one concludes that the age of the universe at T pa 10 10 K would be 
some four years [rather than a few seconds as in standard cosmology]. The 
universe would take some 10 3 years to cool to 10 7 K. With such time periods 
being large in comparison to the free neutron life time, one would hardly 
expect any neutrons to survive at temperatures relevant for nucleosynthesis. 
However, with such a low rate of expansion, weak interactions remain in 
equilibrium for temperatures as low as ~ 10 8 K. The neutron - proton ratio 
keeps falling as n/p ~ exp[—15/T 9 ]. Here T 9 is the temperature in units 
of 10 9 K and the factor of 15 comes from the n-p mass difference in these 
units. There would again hardly be any neutrons left if nucleosynthesis 
were to commence at (say) Jg « 1. However, as weak interactions are still 
in equilibrium, once nucleosynthesis commences, inverse beta decay would 
replenish neutrons by converting protons into neutrons and pumping them 
into the nucleosynthesis channel. With beta decay in equilibrium, the baryon 
entropy ratio determines a low enough nucleosynthesis rate that can remove 
neutrons out of the equilibrium buffer at a rate smaller than the relaxation 
time of the buffer. This ensures that neutron value remains unchanged as 
heavier nuclei build up. It turns out that for baryon entropy ratio i] pa 
5 x 10~ 9 , there would just be enough neutrons produced, after nucleosynthesis 
commences, to give ~ 23.9% Helium and metallicity some 10 8 times the 
metallicity produced in the early universe in the standard scenario. This 
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metallicity is of the same order of magnitude as seen in lowest metallicity 
objects. 

The only problem that one has to contend with is the significantly low 
yields of deuterium in such a cosmology. Though deuterium can be produced 
by spallation processes later in the history of the universe, it is difficult to 
produce the right amount without a simultaneous over production of Lithium 
[19]. However, as pointed out in [1], the amount of Helium produced is quite 
sensitive to i] in such models. In an inhomogeneous universe, therefore, one 
can have the helium to hydrogen ratio to have a large primordial variation. 
Deuterium can be produced by a spallation process much later in the history 
of the universe. If one considers spallation of a helium deficient cloud onto 
a helium rich cloud, it is easy to produce deuterium as demonstrated by 
Epstein [19] - without overproduction of Lithium. 

Interestingly, the baryon entropy ratio required for the right amount of 
helium corresponds to £l b « 0.2. Here il b is the ratio of the baryon density 
to a "density parameter" determined by the Hubble constant: Qt, = p b / p c = 
SnGph/SH^. « 0.2 closes dynamic mass estimates of large galaxies and 
clusters [see eg [20, 21]]. In standard cosmology this closure is sought to be 
achieved by taking recourse to non-baryonic cold dark matter. Thus in a 
linearly scaling cosmology, there would be no need of non-baryonic cold dark 
matter to account for large scale galactic flows. 

b) Density Perturbations 

We conjecture that the universe starts with a homogeneous isotropic dis- 
tribution (eg. a Milne ball) and that an appropriate averaging procedure that 
correctly takes care of retarded effects in a relativistic gravity theory that 
accounts for a creation event, would support the large scale linear coasting 
with the metric described by eqn(l) with K = —1. Perturbations around 
this background shall be assumed to described by 

-SvGST" = 5G, V (3) 

with 5G^ U determined by perturbations around eqn(l). It would be fair to call 
this conjecture as a "Newtonian Cosmology conjecture" as Newton felt that 
one could justify the cancellation of gravity in a homogeneous and isotropic 
universe. Bertschinger [22] and Gibbons [23] have shown how consisent cos- 
mological dynamics, indistinguishable from dynamics in the standard model, 
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follows from a Newtonian proposal in comoving rather than Minkowskian co- 
ordinates. In effect this amounts to stating that a Homogeneous and isotropic 
universe with an arbitrary equation of state can be described by eqn(l), with 
a(t) —t, K = — 1 and that perturbations around it are determined by eqn(3). 
This further amounts to rejecting the unjustified assumptions of continuum 
limit (fluid approximation) and averaging that are usually taken for granted. 
Standard Einstein equations may well hold at the microscopic level but not 
as global averaged equations in a universe with a definite creation event. 
Further, as stated in the introduction, there are alternative gravity models 
that justify the above conjecture. In this article we have restricted ourselves 
to the consequenses of such a conjecture. 

To explore whether such small perturbations in the above model can 
evolve to give structures at large scales, we express the metric as 

ds 2 = {o) g^dx^dx v + 5g^ u dx fM dx v (4) 

with the background metric ^g^ expressed in comoving coordinates as: 

^g^dx^dx" = dt 2 - a(t) 2 ^ ij dx i dx j = a 2 (r])(drj 2 - ^y ij dx i dx j ) (5) 

i] is the conformal time (drj = a~ x dt, a = t = toe 11 ), and 

ll3 = 5 l3 [l-\{x 2 + y 2 + z 2 )}- 2 (6) 

The 5g^ u describe the perturbation. These can be decomposed into scalar, 
vector or tensor type, depending on the way they transform on a constant 
r) hypersurface. In the following, we consider scalar perturbations only. 
Following standard arguments, vector perturbations can be shown to de- 
cay kinematically in an expanding universe and tensor perturbations lead to 
gravitational waves that do not couple to energy density and pressure inho- 
mogeneities [24]. Any scalar perturbation can be decomposed in terms of 
eigenmodes of the laplacian on the constant rj surface: 

V 2 Q = = -k 2 Y (7) 

where '|' represents a covariant derivative with respect to the three metric 7^. 
The most general form for the line element with scalar metric perturbations 
in the longitudinal gauge is (see eg. [25, 26]) 

ds 2 = a 2 (r])[(l + 2^)drf - (1 - m^dx^^)} (8) 
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Similarly, perturbations of the stress tensor can be described in terms of 
velocity and density perturbations in the longitudinal gauge [v^ lon9 \ 5( /on #)). 
We define the matter gauge invariant variables as usual [26] (the dot "'" is a 
derivative with respect to rf): 

V = v {lon9) (9) 
D g = 5 {lon9) + 3(1 + «;)$, w = p/p (10) 
D = 5 ilon9) + 3(l + w )~ (11) 

(X K 

For a perfect fluid, $ = — ^ and the full compliment of scalar perturbation 
equations read: 

AnGa 2 pD = (k 2 + 3)$ (12) 
4irGa 2 (p + p)V = £;((-)*-$) (13) 

(X 

Dg + 3{c 2 s - w)D g a/a + (1 + w)kV + 3wTa/a = (14) 

V + (1 - 3c 2 )V^ = fc(* - 3c^) + ^-£> fl + (15) 
s a s 1+w 9 1 + w 

where c s is the sound speed and r = II; — c 2 5/w vanishes for adiabatic 
perturbations (see eg. [26, 25]). For pure dust w = c 2 = p = and the above 
equations reduce to: 

D g + W = (16) 
V + V- = kV (17) 

(X 

-k 2 ^ = 4ivGa 2 p(D g + 3(tf + (-)^)) (18) 
In terms of the present density, 

C = 47tG Po ^ = ^^p = ^ 6 » 0.3 for Sl b » 0.2 
z oii„ z 

the density perturbation equation simply reduces to: 

\{k 2 + 3)g + 3] + [(A; 2 + 3)g + 2}D g - D g k 2 = (19) 
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k — 1 corresponds to the Hubble scale which is the same as the curvature scale 
in this model. At a redshift ~ 1000, this scale subtends an angle roughly 
.25 degrees. Using constraints from microwave background anisotropy at 
these angles gives D g rs 10 -5 at these scales at the last scattering surface. 
It is easy to see that modes k ~ 1 do not grow. At smaller angular scales 
(large k), the observed anisotropy is expected to fall to much lower values 
[27]. Diffusion damping dampens anisotropies at angular scales smaller than 
about one minute. However, for such large values of k, D g has rapidly growing 
solutions. Eqn.(19) goes as 

D g + D g - Ce^Dg = (20) 

This has exact solutions in terms of modified first and second type bessel 
functions I±, K\. 

D g = Ci(Ce-")3/ 1 ((4Ce-'')3) + ^(Ce'^ ^((ACe'^) (21) 

For large arguments, these functions have their asymptotic forms: 

h — > — 1 exp[2(Ce^n K x — » l ° 6 ± exp[-2(Ce^} (22) 

Even if diffusion damping were to reduce the baryon density contrast to 
values as low as some 10~ 15 , a straight forward numerical integration of 
eqn(19) demonstrates that for k > 3000 the density contrast becomes non 
linear around redshift of the order 50. 

In contrast to the above, in the radiation dominated epoch w — (? s — 1/3. 
In adiabatic approximation, eqns(12,15) imply: 

r . , o ^ \ 3 3C n " 3C • r k ~\~ 3 C t . . 

" + 'iF + 2kh^ " + We^ " + '~ 8 2^ " < 23 > 
For 7] large and negative, small k pertubation equation reduces to 

3D g + QD g - k 2 D g = (24) 

It is straightforward to demonstrate that perturbations bounded for large 
negative 77, damp out for small k. It is more instructive to work with the 
equation for the potential for large k perturbations: 

$ + 3(1 + C 2)$ + [2(1 + 3c 2 s ) + k 2 c 2 s }$ = (25) 
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when 

2e 2r > ■ 

§ long = [3$ + k 2 $] (26) 

For c s = 1/3, the solutions are 

I / ikrj \ 

$ oc e- 2r >e H ~^ (27) 

5 lon9 oc (A; 2 ± s/Zik)e ±{ ^ ] (28) 

Thus we get oscillating solutions for the potential as well as the density 
fluctuations. These fluctuations shall manifest themselves as "acoustic" fluc- 
tuations in the Cosmic microwave background sky to which we shall return 
later. 

We summarize here by noting that fluctuations do not grow in the radia- 
tion dominated era, small k (large scale) fluctuations do not grow in the mat- 
ter dominated well, however, even tiny residual fluctuations O(10 15 ) 
at the last scattering surface for large values of k > 3000 in the matter dom- 
inated era, can in principle grow to the non linear regime. Such a growth 
would be a necessary condition for structure formation and is not satisfied in 
the standard model. In the standard model, cold dark matter is absolutely 
essential for growth. Thus a linear coasting model does not need such dark 
matter to have small density perturbations approach a non linear regime 
necessary for a structure formation theory. 

c) The recombination epoch 

Salient features of the plasma era in a linear coasting cosmology have 
been described in [28, 27]. Here we reproduce some of the peculiarities at the 
recombination epoch. These are deduced by making a simplifying assumption 
of thermodynamic equilibrium just before recombination. 

A recombination process that directly produces a Hydrogen atom in the 
ground state releases a photon with energy B = 13.6eV in each recombi- 
nation. riy(B), the number density of photons in the background radiation 
with energy B, is given by [see eg. [29]]: 

nJB) 16tt - ,-B, 2 x 10 6 ,-13.6, 

= T 3 exp(— ) » ——ex P { ) (29) 

71/ TXi -Lit T 
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Where r is the temperature in units of eV. This ratio is unity at r « .9 
for Q B h 2 ~ 1 and decreases rapidly at lower temperatures. Any 13.6 eV 
photons released due to recombination have a high probability of ionizing 
neutral atoms formed a little earlier. This process is therefore not very 
effective for producing a net number of neutral atoms. The way out of this 
empasse is the two photon emmission. The 2S — > IS transition is strictly 
forbidden at first order. The conservation of both the energy and angular 
momentum in the transition can take place only by emmission of a pair of 
photons. This gives a mechanism for transferring the ionization energy into 
photons with A > \L ya - Since the reverse process does not occur at the same 
rate, this is non-equilibrium recombination. Being a process of second order 
in perturbation theory, this is slow (lifetime ~ 0.1 sec). As recombination 
has to pass through this bottleneck, it proceeds with a rate different from 
the Saha preediction. 

In the red-shift range 800 < z < 1200, the approximate fractional ioniza- 
tion in a linearly coasting cosmology can be approximated by [27]: 

*■ = ^T^^o' 12 ' 25 < 30 > 

The optical depth for Thompson scattering is then: 

dt 



r 7 



ft rz flf 

/ n B (t)x e (t)a T dt = - n B (z)x e (z)a T (—)dz (31) 
jo jo az 

With n b (z) = f]n^{z) = rj x 421.8(1 + zf cm~ 3 , and 

dt 1 



dz~ H (l + z) 2 



(32) 



one can find the red-shift at which the optical depth goes to unity. With the 
residual ionization x etres , we get 

^°- 55 W"' 25 < 33 > 

From this optical depth, we can compute the probability that a photon was 
last scattered in the interval (z, z + dz). This is given by: 

P(z) = e-A „ 7.85 x 10- 3 ( T7 ^) 13 - 25 e^[-0.55( T ^) 14 - 25 ] (34) 
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This P(z) is sharply peaked and well fitted by a gaussian of mean redshift 
z ~ 1037 and standard deviation in redshift Az ~ 67.88. Thus in a linearly 
coasting cosmology, the last scattering surface locates at redshift z* = 1037 
with thickness Az ps 68. Corresponding values in standard cosmology are 
z = 1065 and Az w 80 [30]. 

An important scale that determines the nature of CMBR anisotropy is 
the Hubble scale which is the same as the curvature scale for linear coasting. 
The angle subtended today, by the Hubble radius at z* = 1037, is determined 
by the angular diameter - redshift relation: 

. e . A9)(i + zr) 

r R sm- = smh[ — ] (35) 

For d{9) = 2d H {t R ) = 2H{t R )' 1 = 2[H (1 + z R )}- 1 and r R w (1 + z)/2, this 
gives: 

(I±£«)| = sinh(l) (36) 

or $h ~ 15.5 minutes. 

In standard cosmology, the sound horizon is of the same order as the 
Hubble length at recombination. The Hubble length determines the scale 
over which physical processes can occur coherently. Thus, as we suggest in the 
next section, one expects all acoustic signals to be contained within an angle 
of the order of the angle subtended by the Hubble length at recombination. 

In a linear coasting, the Hubble length is precisely the inverse of the 
curvature scale. However, the sound horizon (s*) is much larger. Strictly 
speaking, the particle as well as the sound horizon are infinite for a linear 
coasting cosmology. As we would see, the CMB anisotropy pattern frozen 
on the Last scattering surface depends on the distance travelled by pressure 
waves from the epoch baryon density starts becoming comparable to radia- 
tion density. The photon diffusion scale is determined by the thickness of the 
LSS. With z* 1037 and Az ps 68, the photon diffusion scale projected on 
the LSS corresponds to an angular size roughly one thirteenth of the Hub- 
ble length at the LSS. This subtends an angle of roughly one minute at the 
current epoch. 

The above scales in principle determine the nature of CMB anisotropy. 
The CMB effectively ceases to scatter when the optical depth to the present 
drops to unity. After last scattering, the photons effectively free stream. On 
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the LSS, the photon distribution may be locally isotropic while still pos- 
sessing inhomogeneities i.e. hot and cold spots, which will be observed as 
anisotropies in the sky today (see eg. [31, 32]). 

d) The CMB temperature anisotropy 

The fractional temperature fluctuation 9(t,x, 7) for a black body is de- 
fined as 

46 = ^ (37) 

Pi 

where p 7 oc T 4 is the spatially and directionally averaged energy density of 
the photons. 9 satisfies the collisionless Boltzmann equation [31, 27]: 

^[9 + $](r ? ,x, 7 ) = $ + ^ (38) 

Before recombination photons and baryons were strongly coupled via 
Compton scattering. In this case the collisional Boltzmann equation is given 
by 

^(9+$) = e+fA(e+$)+fAe-* = f(9 -9- 7 ^+l 7i7j n?) 

(39) 

where f = x e n e ax is the differential optical depth to Thomson scattering 
with <jt = 8na 2 /3ml as the Thomson scattering Cross-section, x e and n e 
are the ionization fraction and total electron density respectively. 9o = 5 7 / 4 
is the isotropic component of 9, 7* = x\ and the quantities 11^ are the 
quadrupole moments of the energy distribution given by 

n? = / ^( W - 7^)49(77, x, 7) (40) 

Perturbations are expanded in terms of eigenfunctions of the Laplacian. 
To effect this, spherical coordinates are defined by the the line element 

dl 2 = lij dx i dx j = -K~ l [d X 2 + smh 2 X (d6 2 + sm 2 6d(f) 2 )} (41) 

where the distance is scaled to the curvature radius x — \J —Kr\. The lapla- 
cian of any scalar function Q is given by: 

7 *Q, y = -Xsinh- 2 X [A (sinh2x ^) + sin^^sin^) + sin" 2 ^] 

(42) 
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Since the angular part is independent of curvature, we may separate variables 
such that Q = XKxWiO, 0) and v 2 = -{k 2 /K + 1). The spherically 
symmetric I = function turns out to be: 



The higher modes are explicitely given by [33] 



(44) 



which reduce to ji{krj) in the flat space limit, where 

Mi = U\ I=0 K V 
K = l 

K x = 1 - (I 2 - l)K/k 2 , I > I. (45) 
All these reduce to unity as K — > 0. The functions are normalized as [31, 34] 

J XU X )XUx)sinh 2 X d X = ^5(u - v')6(l - I') (46) 

These functions can also be generated from their recursion relations. One 
particularly useful relation, which is used in the collisionless Boltzmann equa- 
tion, is [33] 

= -J—kK^Xir 1 + ±±LkKlllxj+ 1 (47) 
dr] v 2/ + 1 ' " 21 + 1 1+1 v K J 

Vectors and tensors needed in the description of the velocity and stress 
perturbations can be constructed from the the covariant derivatives of the 
scalar function Q and the metric tensor, 

Qi = k Q\i 

Qij = k- 2 Q [tJ + -^Q (48) 

where the indices are to be raised and lowered by the three metric 7 U and 7^. 
These vectors and tensors can be used to represent dipoles and quadrupoles, 
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G\ = YQi and G2 = \l l l^Qij- In general any function of position x and an- 
gular direction 7, e.g. the radiation distribution function, can be represented 
as [35] 

00 

%7)=EE*(M,k) (49) 
k 1=0 

where 

G,(x, 7 ,k) = (-A;)- i Q| n ... ii (x,k)P; i - i; (x, 7 ) (50) 
and k 2 = K(k 2 /K + 1) and 

Po = 1, P{ = l l 

^ = ^(37V-7 ij ) 

P»l-'l _ ^ p'Z-'i+l) _J__,(U12 p*3-«i+l) /r-i\ 

^+1 ~ 7TT 7 ^ "ITT 7 ( j 

with parentheses denoting symmetrization about the indices. 
The multipole decomposition of 6 reads: 

8(77, x, 7) = £)e,(i7)Mr 1/2 G,(x,7) (52) 
1=0 

Recursion relations for Gi give the standard hierarchy of coupled equations 
for the /—modes: 

0i = k[Q + $ - W 2 /2 Q 2 ] - f(0 1 - V b ), 
5 

e 2 = k[^K 1 2 /2 Q 1 - 3 -K^Q 3 ] - Af9 2 , 

©/ = fc [^rr^ 1/2 ^-i - 4rl^ 0i+i] " f Qh {l > 2) (53) 

where 7iuj(x) = F 6 Gi(x, 7 ) and if, = 1 - (I 2 - l)K/k 2 . 

Before recombination, Compton scattering transfers momentum between 
the photons and baryons. From the conservation of momentum of the photon 
baryon fluid, we get 

5 b = -kV b + 34/ 
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V b = --V b + k$ + f(V r -V b )/R (54) 

(X 

where R = 3p b /4p 7 . The baryon continuity equation can also be combined 
with the photon continuity equation to obtain [31] 

S b = -k(V b - V 7 ) + ^ 7 (55) 

The differential optical depth f is high enough before recombination mak- 
ing Compton scattering extremely rapid and effective. Since photon-baryon 
tight coupling approximation holds, the Boltzmann eqn(53) and the Euler 
eqn(54) for baryons can be expanded in the Compton scattering time f -1 . 
To zeroth order, we get the tight-coupling identities, 

e = \k 

6i = V 1 = V b 

&i = I > 2 (56) 

These equations merely express the fact that the radiation is isotropic in the 
baryon rest frame and the density fluctuations in the photons grow adiabat- 
ically with the baryons. Substituting the zeroth order solutions back into 
equations (53) and (54), we obtain the iterative first order equations: 

e„ = - j8i + * 
e ^-TTR e ^TTR ke ° + kt (57) 

where we have used the relation R = (a/a)R. The tight coupling approxima- 
tion eliminates the multiple time scales and the infinite hierarchy of coupled 
equations of the full problem. These equations can be combined to form a 
single second order equation, 

®o + + k2 ° 2 s e ° = F ( 58 ) 

1 + JrC 



16 



where the photon-baryon sound speed is 



2 _^7_ = (5g) 

s p^ + pb 31 + R K ' 



assuming that pb ~ and 

f =* + TT5*-T* (60) 

is the forcing function. 

Essential features of CMB anisotropy follow from the above analysis. In 
the rest of this article we shall be content with making rough estimates 
of CMB anisotropy peak locations to judge whether there is any a-priori 
discordance with observations. 

Let us first consider ignoring the time dependence of the potentials \l/ and 
$ and also the baryon- photon momentum ratio R. Then eqn(58) reduces 
to: 

O + k 2 c 2 s & = (61) 



3 

This is a simple harmonic oscillator under the constant acceleration provided 
by gravitational infall and its solution is: 

<3>o(v) = [©o(0) + (1 + J R)$]cos(kr s ) + -!e (0)sin(kr s ) - (1 + R)$ (62) 

KC S 

where r s is the sound horizon given as 

v 

c s drj' (63) 

vo 

where 770 is the epoch of birth of pressure waves. We may take this to be the 
QGP phase transition epoch Tqgp ~ 10 12 K, or even earlier. The following 
results are independent of this initial epoch (cut-off). The distance a sound 
wave travels from the initial epoch till the epoch rj is r s . O (0) and O (O) 
govern the form of the acoustic oscillations and govern the adiabatic and 
isocurvature modes respectively. Eqn(62) also implies, through the photon 
continuity eqn(57), that 

Qi(v) = 3[©o(0) + (1 + i?)$]c s sin(kr s ) + 3ie (0)cos(kr s ) (64) 
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In eqns(62, 64) lie the main acoustic and redshift effects which dominate 
primary anisotropy formation. In the early universe, photons dominate the 
fluid and R — > 0. In this limit, the oscillation takes on an even simpler form. 
For the adiabatic mode, 6 (0) = and Q (v) = [@o(0) + $]cos(kr s ) - $. 
This represents an oscillator with a zero point which has been displaced by 
gravity. The zero point represents the state at which gravity and pressure are 
balanced. The displacement — $ > yields hotter photons in the potential 
well since gravitational infall not only increases the number density of the 
photons but also their energy through gravitational blueshifts. 

However, photons also suffer a gravitational redshift from climbing out 
of the potential well after last scattering. This precisely cancels the — $ 
blueshift. Thus the effective temperature perturbation is Qo(i]) + $ = 
[@o(0) + $]cos(kr s ). This takes care of two sources of primary anisotropy: 
gravity and intrinsic temperature variations on the last scattering surface. 
The final source in this (tight coupling) approximation is the Doppler shift. 
The bulk velocity of the fluid along the line of sight is V 1 j \[Z = Qi/y/3. This 
causes the observed temperature to be Doppler shifted. Eqn(64) shows that 
the acoustic velocity Gi is tt/2 out of phase with the temperature. Because of 
its phase relation, the velocity contribution fills in the zeros of the tempera- 
ture oscillations. The rms temperature fluctuation is identically smoothened 
out leaving no residual primary anisotropies. 

This null result gets altered once we take into account of a non-vanishing 
R. Baryons add gravitational and inertial mass to the fluid, m e ff = 1 + R, 
without affecting the pressure. This decreases the sound speed and changes 
the balance of pressure and gravity. Gravitational infall now leads to greater 
compression of the fluid in a potential well, i.e. a further displacement of 
the oscillation zero point. Since the redshift is not affected by the baryon 
content, this relative shift remains after last scattering to enhance all peaks 
from compression over those from rarefaction. To demonstrate the effect of 
R 0, let us assume it to be constant and non-vanishing after an epoch 
r] = rji and negligible earlier. For an adiabatic perturbation, 6 (0) = we 



have: 




(65) 




k 



k 



(66) 
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Thus after 77 = 771 , the solution reads: 

8(77) + (1 + R)$ = [6(0) + $][cos(-^Ui)cos( k (77 - 771)) 

V3 ^/3(1 + R) 



k k k 

1 + fl)sin(-r=77i)sin( =(77 - 771))] + R$cos(^^==(r/ - 771)) 

V3 ^3(1 + R) 73(1 + R) 



(67) 

The phase of the oscillation freezes at the last scattering surface 77 = 
77* and determines the observed fluctuation pattern. For a linear coasting 
cosmology, 1 771 1 is large. It is in fact infinite for a universe starting at t — 0, 
and can be made finite (but large) only by choosing a lower cutoff for t. The 
angular power spectrum of the above expression is obtained by transforming 
it with respect to the ultraspherical bessel functions [31, 34]. A given k mode 
would contribute to two distinct angular scales: (1) A very small angular scale 
corresponding to the large phase of the first two harmonic functions, and (2) 
the angular scale corresponding to the third harmonic function: 

+.R$cos(-.= (77* - 771)) 

>/3(l + R) 

For PS 0.2 and h ~ 0.7, significant deviation of c s from l/y/3 occur around 
the matter radiation equality at redshift given by l + z eq ps 3.9 x 10 4 (Q b h 2 ) ~ 
4000. Thus only those fluctuations that originate at redshift z = few x z eq 
would give rise to observable temperature anisotropies at LSS which locates 
at z* ~ 1037. The argument of the above harmonic pattern freezes at LSS 
and would give rise to "peaks" at 

a Rec \n(— ) = n7r (68) 



The corresponding wavelengths are at 



7T 



K = Y ( 69 ) 



The first "peak" is at k a = and the second at roughly 9 ps 15 minutes. 
Compression peaks occur at odd values of n and are located at angles 9% = 
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15/ 'mr minutes. The entire analysis can be repeated for isocurvature modes 
to find the peaks located at even n at l n — 15/(n + \) minutes. 

The above analysis has not taken into account the variation of R. For 
a qualitative comparison one has to integrate eqn(58). Using a WKB ap- 
proximation (which is good for modes for which the time scale of variation 
of sound speed is greater than the oscillation period) it can be seen [31] that 
the phase structure only suffers a minor change from the above analysis while 
most of the effect of changing R effects the amplitude of the solution. 

The exact profile of the anisotropy would be determined by the choice of 
the nature of initial conditions (adiabatic or isocurvature), the chosen initial 
power spectrum, and the growth of perturbations after z* (these determine 
the late or the integrated SW effect), aspects of reionization etc. As photons 
possess a mean free path in the Baryons A c ~ f _1 due to Compton scattering, 
the photon-baryon tight coupling breaks down at the photon diffusion scale. 
Finally, after last scattering, photons free stream toward the observer on ra- 
dial null geodesies and suffer only the gravitational interactions of redshift 
and dilation. Spatial fluctuations on the last scattering surface are observed 
as anisotropies in the sky. Free streaming thus transfers I = inhomo- 
geneities and I = 1 bulk velocities to high multipoles as the /—mode coupling 
of the Boltzmann eqn(53) suggests. Isotropic / = density perturbations are 
thus averaged away collisionlessly. 



3 Summary 



The main point we make in this article is that in spite of a significantly 
different evolution, the recombination history of a linearly coasting cosmology 
can be expected to give the location of the primary acoustic peaks in the same 
range of angles as that given in Standard Cosmology. Given that none of 
the alternative anisotropy formation scenarios provide a compelling ab initio 
model [39] , it is perhaps best to keep an open mind to all possibilities. As the 
large scale structure and CMB anisotropy data continue to accumulate, one 
could explore the general principles for an open coasting cosmology to aid in 
the empirical reconstruction of a consistent model for structure formation. 
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Finally, we are tempted to mention that a linear coasting cosmology 
presents itself as a falsifiable model. It is encouraging to observe its con- 
cordance !! In standard cosmology, falsifiability has taken on a backstage 
- one just constrains the values of cosmological parameters subjecting the 
data to Bayesian statistics. Ideally, one would have been very content with 
a cosmology based on physics that we have already tested in the laboratory. 
Clearly, standard cosmology does not pass that test. One needs a mixture of 
hot and cold dark matter, together with (now) some form of dark energy to 
act as a cosmological constant, to find any concordance with observations. 
In other words, one uses observations to parametrize theory in Standard 
Cosmology. 

Linear coasting presents a very distinguishable cosmology. Recombina- 
tion occurs at age ~ 10 7 years as opposed to ~ 10 5 years in standard cosmol- 
ogy. The Hubble scale at decoupling is thus two orders of magnitude greater 
in linear coasting. This fact, coupled with the absence of any horizon, could 
well have falsified linear coasting. Any concordance with observations is 
therefore very significant. The model comes with its characteristic predic- 
tions. Thus linear coasting has the potential of relegating the need for any 
form of dark matter or dark energy (or for that matter, any physics not al- 
ready tested in the laboratory) to the physics archives where they enjoy the 
same status as ether and phlogiston. The message this article is to convey is 
that a universe that is born and evolves as a curvature dominated model has 
a tremendous concordance and there are sufficient grounds to explore models 
that support such a coasting. 
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Figure 1: Hubble diagram, taken from the supernova cosmology project. 
"The curve for (fl^, Q M ) = (0, 0) is practically identical to the best fit un- 
constrained cosmology" [13]. 
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